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The problem of two interacting particles moving in a d- 
dimensional billiard is considered here. A suitable coordinate 
transformation leads to the problem of a particle in an un- 
conventional hyperbilliard. A dynamical map can be readily 
constructed for this general system, which greatly simplifies 
calculations. As a particular example, we consider two iden- 
tical particles interacting through a screened Coulomb poten- 
tial in a one-dimensional billiard. We find that the screening 
plays an important role in the dynamical behavior of the sys- 
tem and only in the limit of vanishing screening length can 
the particles be considered as bouncing balls. For more gen- 
eral screening and energy values, the system presents strong 
non-integrability with resonant islands of stability. 



A system of two interacting bodies moving in an oth- 
erwise free space, is one of the few integrable problems 
known. The reduction to the one-body central force 
problem allows a solution by quadratures However, 
once the translational symmetries are broken, as when 
the system is placed inside a billiard, the center-of-mass 
(CM) and angular momenta are in general no longer con- 
stants of motion. In this case, the classical dynamics of 
the system may be chaotic even when the geometry yields 
an otherwise fully integrable one-particle case, as we shall 
see below. 

On the other hand, recent experimental realizations of 
billiards, such as suitably shaped resonators and quan- 
tum dots |2|,|3|], have allowed the study of the quantum 
manifestations of well-known classical non-integrability 
in some billiards pL§. In the case of quantum dots, 
disagreement between theory and experiment has been 
attributed to geometrical factors ||. A considerable 
amount of theoretical work exists on the effect that ge- 
ometry has on the integrability of dynamics in billiards 
|,|-|, as well as on their quantum analogs How- 
ever, the possibility of more than one particle in the 
quantum dot leaves the usual one-particle approach in- 
complete. In fact, some experiments have pointed out 
the importance of electron-electron interaction on various 
features observed in such mesoscopic systems [|l0| . In this 
article, we explore the role of the electrostatic interaction 
introduced when two particles are in the billiard. A for- 
malism for billiards in any dimensions is developed, and 
as an example, we apply it to the one-dimensional case. 
Since we are interested in the role of the electrostatic 
interaction in mesoscopic systems, we consider particles 



interacting through a screened Coulomb potential. 

The hyperbilliard. The problem of two point masses 
moving along a finite line and suffering elastic impacts 
with the end walls and between themselves, can be trans- 
formed to the motion of one 'particle' moving in a tri- 
angular billiard. The coordinates of the particle in this 
billiard are the coordinates of the original masses. The 
ratio of the masses determines the integrability of the 
system [fjj, being regularizable for a particle mass ratio 
of 1 and 3 (or ±) §. 

We now introduce an interaction between the particles 
and consider the ci-dimensional case. Let q^, Pi, (i = 1,2) 
be the position and linear momentum of the i-th parti- 
cle. The motion takes place in a <i-dimensional billiard, a 
compact simply-connected region of R d whose boundary 
is denoted by T. We assume that T is piecewise smooth 
and defined by v surfaces, Tj = {q : fj(q_, ctj) — 0},j — 
1, u, where fj and ctj denote the function and the set 
of constants which characterize the j-th surface. These 
functions define subspaces of dimension d— 1 in R d . To fix 
ideas, we restrict ourselves to flat surfaces, i.e., for d = 2 
(3) the billiards are simple polygons (polyhedrons). 

The formalism developed here can be applied to any 
central-force interaction between the particles. We have 
selected the screened Coulomb potential, i.e., the Yukawa 
potential given by ^(qi,q2) = e~ A l q2 ~ qi l/|q2 — qi|, 
where A -1 is the screening length. Notice that this poten- 
tial goes to a S— function when A — > oo. In this limit, the 
particles behave as bouncing hard-core balls, i.e., non- 
interacting impenetrable point particles, for which the 
dynamics can be integrable, as described above. Hence, 
for a given energy, A plays the role of the perturbation 
parameter. Due to the interaction, a finite value of A 
determines the finite effective radius of the particles for 
a given total energy, as described below. 

Considering for simplicity identical-mass particles 
(mi = rri2 = 1), the Hamiltonian for the system is writ- 
ten as 

H = f + f + V( qi , q 2 ) + £ £ Ufofaaj)], (1) 

i=l 2=1 

where the function U[fj(cii,aj)] represents the infinite 
repulsion potential exerted by the j-th hard wall on the 
i-th particle. Analytically, this function could be written 
in terms of Heaviside functions with a large prefactor. 
In practice, the normal component of the velocity of the 
incident particle will be reversed at the moment of bounc- 
ing on the billiard walls. 
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The Hamilton equations can be written as = Pi, and 

Pi = -V g4 V(qi,q2) + Ej=i A j(PiWj{<li,(Xj)], where 
i = 1,2 and the vector function Aj represents the change 
of momentum due to the bounce on the j-th wall. Given 
d, the number of degrees of freedom is 2d. Hence, the 
phase space of the system is 4d-dimcnsional. 

We now introduce a transformation to center-of-mass 
and relative coordinates R = (qi + c^j/M, and r = 
q2 — qi, respectively, where the total mass M = 2, and 
the reduced mass p = |. These equations define a new 
space of coordinates p = (r, R), which is 2d-dimensional. 
In this space, we have a new set of equations for the 
boundary of the billiard, say Fj(p, ocj), j = l,...,v. Every 
function Fj now defines a subspace of 2d — 1 dimensions 
in p-space. 

The Hamilton equations are transformed then to r = 
p/fi, R = P/M, and 

p = -V r V(r) + A j (P- p ) ^( r > R a i)\ > 

3=1 



P = ^B,(p,P)^(r,R,a,)]. 



(2) 



As before, Aj and Bj represent the change of the mo- 
menta p and P, respectively, due to the bounce on the 
j-th wall. 

Notice that these equations describe the motion of one 
particle in the p hyperspace, i.e., we have constructed the 
hyperbilliard. The description of a system composed by a 
few masses in terms of one particle in a hyperspace has 
been used for several cases, including billiards |^,0,|l2|. 
Usually, the hyperspace is constructed without introduc- 
ing transformations of the coordinates. Here, however, 
the change to the CM coordinates allows one to get a 
map in a simple way. 

Notice that bounces of the particle in the hyperbilliard 
correspond to bounces of the masses in the real billiard. 
The walls of the billiard cause the breaking of the trans- 
lational symmetry of the system, and as a consequence, 
the CM momentum is no longer a constant of motion. In 
the case of non-interacting and equal-mass particles, the 
changes in the CM momentum are determined only by 
the geometry of the billiard. In our case, however, the 
interaction couples the CM and relative momenta after 
a bounce, which in turn depend on the momenta of each 
of the original masses. The rotational symmetry is also 
broken in general and the generator of rotations is no 
longer a constant of motion either. 

The map. Hamilton equations in p-space indicate that 
between bounces the particle moves freely along the CM 
coordinate whereas the central force V(r) acts only along 
r. The motions are independent, and only become cor- 
related at each bounce, as the corresponding momenta 
are changed while keeping the total energy constant. We 
take advantage of this fact: Consider that the particle at 
the n-th bounce has the coordinate p n — (r„,R„). The 
condition that the time spent by the particle until the 



next bounce on the j-th wall be the same along the r 
and R coordinates, 



-(pnoPn+l) — TR k {p n , Pn+l), 



k = l,...,d, (3) 



represents an interesting opportunity. Here, r r (rn k ) 
refers to the time along the relative (fc-component of 
CM) coordinate. The times TR k for the free motion be- 
tween collisions can be calculated easily. The l.h.s. in 
(0) can be obtained by noting that the motion along r 
becomes separable and the time t t can then be calcu- 
lated by quadratures, as illustrated below. Equations (||) 
and the equation corresponding to Fj result in a set of 
nonlinear algebraic equations for p n +i- We call this set 
the map of the billiard since it indeed expresses p n +i in 
terms of p n . This procedure can be easily carried out at 
least formally in the general case. Notice that this map 
has not been obtained by means of the usual linearization 
procedure |11|, but rather as an extension of Benettin's 
procedure p|]. The ID case, explained in detail now, 
provides a clear example of this procedure. 

The ID billiard. This system is defined by walls at 
the end points q = ±| . Because of the interparticle 
repulsion, the particle 1 (2) never reaches the boundary 
2 (1). This implies that Aj, here associated with bounces 
on the j-th wall, will describe bounces of the j-th particle 
only. The Hamilton equations for the r — R coordinates 
are then p = pf, P = MR, and 



. _ dV(r) 
^ dr 

P = B l 8{R- 
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According to the arguments of the 8- functions, the point 
boundaries are transformed into lines in p-space, which 
define a billiard with an isosceles-triangle shape, simi- 
lar to the case of non- interacting hard core particles ||, 
although here is in the r — R space. The base of this 
triangle in our case acts as a repulsive wall of potential 
V{r). The closest approach to the repulsive wall by the 
particle (the turning point), depends on the energy asso- 
ciated with the relative motion, e = E — P 2 /2M. 

The functions giving the change of momentum are sim- 
ple. For example, for bounces of the i — th particle (on 
i — th wall) we have Ai = — p ± 2pP/M, where P and p 
are the momenta before the collision, and the + (— ) sign 
refers to i = 1 (2). 

For a pure Coulomb potential (A = 0), r r can be calcu- 
lated analytically, so Eq. (||) can be written in the form 



T{T e (p n ),T e (p n+1 )} = 



Rn+l — Rr, 



P/M 



(5) 



where r is the time elapsed going from p n to p n +i, ex- 
pressed in terms of the time T spent by the particle from 
the turning point to p 



T e (^f+(^cosh-W. 



(6) 
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For A < 1 we can expand V(r) to first order and obtain 
the same expression, except that E is shifted to E—X. For 
all different initial conditions there are only a few possible 
trajectories which can be determined by analyzing the 
momenta. A simple algorithm can then be obtained to 
determine the Poincare surfaces of section. (The details 
of the motion in p-space will be presented elsewhere.) 
This nontrivial algebraic map provides a full description 
of the dynamics. Its use simplifies calculations a great 
deal, and allows one to better characterize the system, 
as we describe below. 

To characterize the dynamics, we determine the 
Poincare section (PS) at a phase such that one of the 
masses is fixed, say, as it just bounces on the wall. Then 
we plot the position and momentum of the other mass. 
Because of the indistinguishability of the particles, the 
topology of the PS does not depend on which mass is 
selected. In fact, the surfaces are identical, except for 
left-right exchange symmetry. 

Between bounces on the walls, the masses approach 
each other a distance given by their relative energy. The 
shortest distance of approach r m , can be considered as 
twice the minimum effective radius of the particles (for 
zero CM energy). This effective radius is a characteristic 
of the system and its dependence on the total energy 
E and the inverse screening length A is shown in Fig. 
[I], as obtained from the condition E = V(r m ). When 
A > 1, the interaction is short-ranged, which results in 
nearly-free particles for some moments. Note that for all 
values of E and A, the initial condition = — i q^ = 

\,p^ = \/[E — V(1)],P2 = — Pi ' corresponds to a 
periodic orbit and we call it the symmetric motion. 

We now fix E — 1.56 and change A. Fig. [2^, shows the 
PS for the "short-range potential" case, A = 20. These 
results were obtained by direct numerical integration of 
the equations of motion. Here q± = — i, which means 
that we are plotting the position and momentum of the 
particle 2 as particle 1 is at its (left) edge of the billiard. 
According to the KAM theorem some invariant tori 
will be preserved under the interaction, although they 
are somewhat deformed. These periodic or quasiperiodic 
orbits lie inside the primary islands of stability, which in 
our case are situated around the fixed point correspond- 
ing to the symmetric motion, the latter one represented 
by means of x . Higher values of A present a similar PS, 
but the chaotic region fills more and more of the avail- 
able space consistent with E fixed. It is clear that only 
for the case of point particles, i.e. zero screening length 
(1/A = 0), the bouncing-ball behavior is observed. It is 
also possible that the infinite energy limit (with r m = 0) 
would be similar, as Fig. 1 also suggests. 

Secondary islands appear for A ~ 1. These islands are 
due to the inter-particle interaction (]ll| and correspond 
to correlated motion, as when for example, particle 1 
bounces twice and the other once. For each pair of val- 
ues (A, E) there is a specific island structure. The orbits 
in the secondary islands become unstable for the "short- 



range" case because for some instants the particles are 
nearly free, the memory of the previous motion is lost, 
and the correlation is destroyed. Hence, as A decreases 
(Fig. |^b), the number of stability islands increases. In 
this case qi = i, so that the PS shows the position and 
momentum of particle 1 when particle 2 is at the right 
edge of the billiard. Notice that the available region of 
space decreases for smaller A, as the total energy is fixed 
and the inter-particle potential energy has a stronger con- 
fining effect. 

The results using the map described before are now 
presented. Figure [| shows the PS for A = obtained 
by solving (a) the differential Hamilton equations in q- 
space, and (b) the algebraic Eq. (||). The graph obtained 
by means of the latter has been reflected about (72 = h 
for easy comparison. The two A = PS are topologically 
identical (if traversed in different sequences) . The agree- 
ment is excellent even for Awl, while the computation 
time is substantially reduced (~ 10 3 times) if the map is 
used. 

Using the map, we have calculated the Lyapunov ex- 
ponent a for A < 1, following the procedure of Ref . [|l3| . 
Figure || shows that, as A increases, the fraction of phase 
space filled by the chaotic sea increases also. This is re- 
flected in the Lyapunov exponent (not shown), which for 
a constant energy [E — 1.56) increases monotonically 
(from 0.34 to 0.59) with the inverse screening length A 
(0 to 1). Increasing energy produces similar curves with 
ever larger values of a. 

A general formalism for two interacting particles in 
a d-dimensional billiard has been presented. The one- 
dimensional case with a screened Coulomb potential was 
shown to exhibit soft chaos. Only in the case of infinite 
screening length (or energy), the particles can be consid- 
ered as bouncing balls. These results suggest that the 
effects of electrostatic interaction between electrons in 
quantum dots, for example, may play a very important 
role in the quantum-classical correspondence and they 
should be considered when these systems are studied. 
The analysis of the quantum mechanical analog of the 
billiard system described in this work is now in progress. 
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FIG. 1. The closest approach between the particles, r m , 
as a function of the inverse screening length A, for different 
energies E. 



FIG. 2. Poincare sections for E = 1.56. a) qi — — |, 
A = 20, b) qi — i, A = 0.6. Crosses (x) indicate the sym- 
metric periodic motion. 



FIG. 3. Poincare sections for the selected energy E — 1.56, 
A=0 and qi = — a) Solving the Hamilton equations and b) 
using the map. 
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This figure "fig2a.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9703098vl 



This figure "fig2b.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9703098vl 



This figure "fig3.gif" is available in "gif" format from: 
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